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Abstract. We derive the expression of the stress tensor for one and two-component lipid membranes with 
density and composition inhomogeneities. We first express the membrane stress tensor as a function of 
the free-energy density by means of the principle of virtual work. We then apply this general result to 
a monolayer model which is shown to be a local version of the area-difference elasticity (ADE) model. 
The resulting stress tensor expression generalizes the one associated with the Helfrich model, and can be 
specialized to obtain the one associated with the ADE model. Our stress tensor directly gives the force 
exchanged through a boundary in a monolayer with density and composition inhomogeneities. Besides, it 
yields the force density, which is also directly obtained in covariant formalism. We apply our results to 
study the forces induced in a membrane by a local perturbation. 



PACS. XX.XX.XX No PACS code given 



1 Introduction 

Biological membranes consist principally of a lipid bilayer. 

. In each monolayer, the molecules are arranged with their 
hydrophobic ends directed inward toward the second mono- 
layer, and their hydrophilic ends directed outward toward 

■ the surrounding water [T][7j • Membranes are fluid and pos- 
sess a weak resistance to bending, so they fluctuate strongly 
under the effect of thermal excitations [3J. In many rele- 

' vant situations, their total area may be treated as a con- 

, stant: this is due to the existence of a practically constant 
area per lipid and to the fact that lipids are virtually in- 
soluble in water [7j . As a consequence, membranes develop 
a lateral tension in response to external forces [4][5]. 

Although the constant area approximation works finely 

, for many fundamental aspects (see, e.g., Refs. 
[TUlfTTlfTT] ) . there are important situations in which it is 
essential to take into account the variability of the lipid 
density in each monolayer. Indeed, an essential parameter 
in the phase diagram of vesicle shapes is the difference be- 
tween the areas of the two monolayers [T4l[Fol[T6l[18] . and 
the most accurate equilibrium description of the elasticity 
of membranes thus involves a monolayer "area difference 
elasticity" (ADE) term in addition to the usual bending 
terms (171118(115] . There are also spectacular instabilities 
that may be triggered by a local expansion of one mono- 
layer, yielding cristae-like invagination |21j . ejection of 

a Author for correspondance. 



tubules [22] or bursting and curling of polymcrsomes |23) . 
Finally, the dynamics of membranes at the sub-micron 
scale is controlled by the intermonolayer friction result- 
ing from local velocity differences between the monolayers 
fluids that must be treated as compressible [241125] . 

In this paper, we derive a fundamental tool that will 
help to rationalize these phenomena: the stress tensor act- 
ing in a membrane monolayer of variable shape, vari- 
able lipid density and variable composition. Understand- 
ing forces in complex membranes with various degrees of 
freedom is crucial to understand their equilibrium shape 
and their dynamics [251177] . The divergence of the stress 
tensor gives the density of clastic forces in the membrane, 
which is the basis of a dynamical description. It is very 
useful in the study of shape instabilities [75]. Moreover, 
the stress tensor provides the forces directly exchanged 
through a boundary. This is a valuable information that 
can be used, e.g., in the calculation of membrane-mediated 
interactions [751130] . or in the study of the spatial distri- 
bution of forces on the edge of an object embedded in the 
membrane. 

The outline of our paper is the following: In Sec. 2, we 
present a formal derivation of our stress tensor from the 
principle of virtual work. In Sec. 3, following Ref. [75], we 
construct a free-energy density for one and two-component 
monolayers that extends that of the ADE model. In Sec. 4, 
we derive the stress tensor associated with this free-energy 
density, and in Sec. 5 we discuss the associated force den- 
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sity. Finally, in Sec. 6, we show how our results enable to 
understand the forces and the dynamics in the case of a 
local perturbation of a membrane by the microinjection of 
a reagent. 



2 Stress tensor formal derivation 

Let us consider one monolayer of the bilayer. We describe 
its shape in cartesian coordinates by the equation z = 
h(x, y) of a surface S parallel to its hydrophobic interface 
with the other monolayer. Such a description, based on 
the height of the membrane with respect to a reference 
plane, is often referred to as the Monge gauge. We do 
not assume yet that the membrane is weakly deformed. 
Let p(x,y) be the projected mass density, i.e., the lipid 
mass per unit area of the reference plane (x,y). In order 
to study the case of a two-component monolayer, let us 
denote by </>(x, y) the local mass fraction of one of the two 
lipid species, say species number 1. The case of a one- 
component monolayer can be obtained by setting <f> = 0. 

Let f(p, 4>, hi, hij) be the projected free-energy density 
of the monolayer (i.e., the free energy per unit area of the 
reference plane). Here and in the following, Latin indices 
represent cither x or y (not z) and hi = hij = didjh, 
etc. Note that we are assuming that the free energy de- 
pends only on the mass density, on the local lipid compo- 
sition, and on the slope and curvature of the monolayer. 
We have thus neglected the gradients of the curvature and 
of the lipid density and composition. The former approx- 
imation has already proved successful [511TB]: the latter is 
justified by the fact that the correlation length of the den- 
sity fluctuations should not be larger than the monolayer 
thickness (far from a critical point). 

Let us consider an infinitesimal cut with length ds sep- 
arating a region A from a region B in the monolayer, and 
let us denote by m the normal to the projected cut di- 
rected toward region A. The projected stress tensor £ 
relates linearly the force df that region A exerts onto re- 
gion B to the vectorial length m ds of the projected cut 
through 

df = E-mds. (1) 



This defines the six components of the projected stress 
tensor: and S z j, where i € {x,y} and j € {x,y} [32] , 
To determine the projected stress tensor, we shall fol- 
low the method presented in Ref. [32], which is based on 
the principle of virtual work. Let us consider a monolayer 
patch standing above a domain Q of the reference plane. 
This patch is supposed to be a closed system with fixed 
total mass of each lipid species. Its free energy reads 



F = I d r f (p,4>,hi,hij) 
In 



(2) 



In order to deal with arbitrarily deformed states of 
the monolayer patch, we assume that in addition to the 
boundary forces (and torques) exerted by the rest of the 
monolayer, the patch is submitted to a surface density 
w(x, y) of external forces, and to individual external forces 



acting on the molecules and deriving from a one-body po- 
tential energy v a (x, y) for the lipid species a £ {1, 2}. The 
former forces control the shape of the patch and the latter 
control the mass density distribution of the lipids within 
the patch. The total potential energy corresponding to 
these latter forces can be written as 



E p = drtivi + dri2V2 = I c 
Jn Jn 



Ml 



.(3) 

where d 2 r = dx dy, and p a denotes the mass of one lipid 
of the species a. Introducing v = V2I p-2 and u = V\lp>\ — 
v il 1 ^2, this potential energy can be rewritten as: 



E p = d 2 r [pv 
J n 



(4) 



At equilibrium, the lipid density p, the composition 4> and 
the shape h of the monolayer are controlled by the external 
actions represented by u(x,y), v(x,y) and w(x,y). 

Let us study a small deformation of the monolayer 
patch at equilibrium: Q — > [2 + SS2, h — > h + Sh, p — > p+Sp 
and 4> —¥ + 5<f>. Each element of the patch, initially at 
position {x,y), undergoes a displacement 5a{x,y), with 
Sa z = Sh + hkSak [32]. The variation of the free energy of 
the monolayer patch during the deformation reads 



SF = 



d 2 r 



df . df , ± df „ 
op o<p ohi 



dhu 



d 2 rf, 



,)<> 



(5) 



We now perform two integrations by parts, and we use the 
relation 



rf 2 , 



ds rrii 5a,i , 



(6) 



where B denotes the boundary of SI. Assuming that the 
translation of the monolayer edges is performed at a fixed 
orientation of its normal, so that Shj = —hj^Sak along the 
boundary |32| . we obtain 



SF = / d 2 r 



[Op 



+ / ds mi < fScii + 



df . df, ± SF ' 
-^=Sp+ ^-5<j)+ —Sh 



5h 

df g df 

dhi J dha 



df df 



' >J i)h, , dhi ' 



df_ 

dha 



Sa z 
Sa k 



(7) 



where 



Sh dh jk ohi 



(8) 

The potential energy variation during the deformation is 



SE n 



d 2 r [(v + u(f>) Sp + up 1 
ds nii Sai [vp + uptfi] . 



(9) 
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The total variation SF + SE p of the energy of the system 
must be balanced by the work SW done by the surface 
force density w and by the boundary forces exerted by the 
rest of the membrane on our patch. Since the translation 
of the monolayer edges is performed at a fixed orientation 
of its normal, the torques produce no work. We may write 



SW 



d r [w k Sa k + w z Sa z 



/ ds rrii [E ki Sa k + S zi Sa z 

JB 



(10) 



As the monolayer patch is considered as a closed system, 
the total mass of each lipid species in the patch is constant 
during our deformation: 



d 2 rp = M and / d l r , 



M 1; 



(11) 



where M and Mi are constants. Let us introduce two con- 
stant Lagrange multipliers A and \i to implement these two 
global constraints. The relation SF + 5E p - SW + XSM + 
fiSMi = must hold for any infinitesimal deformation of 
the monolayer patch. The identification of bulk terms in 
this relation yields 



S£ 
Sh 



w z = — and w k = -h k w z , 



df 

— = -(v + X) - (U + fJL)<f>, 

df 



(12) 
(13) 
(14) 



By identifying the boundary terms and using (|13|) , we ob- 
tain the components of the membrane stress tensor: 



dk-^J— ) 



Oh 



dh k 



df 
dh kj 

d]_ 

dh 3 



'■hi ■> 



df 
dh 



(15) 



(16) 



which generalizes the result of Rcf. [32] to the case where 
there are inhomogeneities in p and <f>. Note that the frac- 
tion of each lipid species does not appear explicitly in 
this result. Therefore, Eqs. (|T51) - p^|) hold both for one- 
component monolayers and for two-component monolay- 
ers. Note also that £ does not depend directly on the 
external actions, which confirms its intrinsic nature. 

Comparing our result with Ref. [32] shows that tak- 
ing into account lipid density variations only changes the 
isotropic term of the stress tensor, which now reads / — 
pdf/dp. This term is reminiscent of minus the pressure 
of a two-dimensional homogeneous fluid binary mixture 
with area A described by a free energy F(T, A, N\, N2) = 
Af(T,p,<j>): 



p 



dF 
dA 



f 



T,Ni,N 2 



df 
dp 



(17) 



While the last expression can be used locally in a non- 
homogeneous fluid mixture, the case of the membrane is 
more complex since its free energy depends on the cur- 
vature. The "surface pressure" in a membrane is some- 
times defined as pdg/dp — g where g is the part of the 
membrane free-energy density / that depends only on p 
[34] or, equivalcntly, what remains of / for a planar mem- 
brane [35] . Interestingly, we find that the isotropic part 
of the membrane stress tensor does not identify to minus 
this surface pressure, since it is the complete, curvature- 
dependent projected free-energy density / that appears in 

/- pdf/dp. " 

The divergence of the stress tensor gives the force per 
unit area, p, exerted by the rest of the monolayer on the 
patch. By direct differentiation, we obtain: 



Pz 

Pi 



SF 



df , df. 



djEij = -jrhi - pdi— + -^rdifi, 



Sh 



dp dtp 



(18) 
(19) 



where we have usjsd dif = (df /dp)dip + (df /d(j))di4> + 
(dffdhj)hij + (df/dhj k )hij k . At equilibrium, we can use 
([T^]) — ([T^]) to express p, which yields 



Pz 

Pi 



-Wz 

- Wa 



pdiV + 



(20) 
(21) 



These relations constitute the balance of surface force den- 
sities for the monolayer at equilibrium. In particular, p 
vanishes at equilibrium when the membrane is submitted 
to no external actions (i.e. w = and u = v = 0). 



3 Monolayer model 

3.1 Free-energy density in terms of local variables 

3.1.1 One-component monolayer 

Let us derive the elastic free energy of a monolayer in 
a bilaycr from basic principles, first for a one-component 
monolayer. We shall recover and extend the model of Ref. [? 
which is actually a local version of the ADE model. 

We assume that the free-energy density / per unit area 
of the monolayer depends only on the mass density p and 
on the local principal curvatures C\ and ci of this mono- 
layer. As in the previous section, the gradients of the cur- 
vature and of the density are neglected in our description. 
Note that, unlike / and p, f and p are the free energy 
and the mass per actual unit area of the monolayer, and 
not per projected unit area. We use for both monolayers 
the density and the principal curvatures defined on the 
same surface S of the bilayer, so that the curvatures are 
common to the two monolayers. 

We will consider the physically relevant regime of cur- 
vature radii much larger than the membrane thickness. 
We will also restrict ourselves to small variations of the 
density around a reference density pq. Note that it can 
be convenient to take po different from the equilibrium 
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density p cq of a plane monolayer with fixed total mass, 
for instance to study a monolayer under tension. Let us 
define 



P~ Po 



Po 

H=(a + c 2 ) e 
K = cic 2 e 2 = (D(e 2 



0(e), 
0(e), 



(22) 

(23) 
(24) 



where e is a small length in the nanometer range that 
allows to define the scaled total curvature H and the 
scaled Gaussian curvature K. Since we typically expect 
10" 4 < \ r \ < 1CT 2 and 1CT 4 < \ae\ < 10~ 2 , it is sensi- 
ble to assume that r and H arc 0(e) while K is C(e 2 ). 
The free-energy density / is a function of these three non- 
dimensional small variables. 

To study small deformations, we write a second-order 
expansion of /: 



f(r, H, K) = a + A 1 H + A 2 {r- H) J 
+ AiK + Ofe 3 ) . 



A, H 2 



(25) 



Three comments are due here. i)We have not included any 
term linear in r in this expansion. Indeed, the total mass 
of the monolayer, i.e. the integral of p = po(l + r), is 
assumed to be constant, so including a term linear in r 
is equivalent to redefining the constant term <jq. ii) The 
freedom associated with the choice of e allows to set the 
coefficient of —rH equal to twice that of r 2 . iii) All the 
coefficients in (|25[l depend on the reference density pq. 
Wc shall come back in the following on the constitutive 
relation oo(po)- 

Defining the constants k, k, cq and k through 



Ai = - 



KCq 

2e ' 



A 2 = \k, 



A* = 



4e 2 



Aa 



2^' 



(26) 



and setting c = c\ + c 2 , we obtain the following expression, 
which generalizes those of Refs. [231133] : 



k 2 KCo K 2 , K 

J =cr o+2 (r-ec) -— c +^ c +2 ClC2 " 



(27) 



Note that all the above terms have the same order of mag- 
nitude. Indeed, since typically k w 10 _1 J/m 2 , k w k ps 
10 _19 J, e w lnm and Cq 1 « 50 nm (see, e.g., Refs. [H 
[TBH37]). we have A 2 « A 3 w A 4 w 100 Ai. The advan- 
tage of the procedure we have employed is that we control 
precisely the order of the expansion. 

In (|2"T|) . e can be interpreted as the distance between 
the surface S where c\ , c 2 and r are defined and the neutral 
surface of the monolayer [3"Tll38j . As a matter of fact, the 
density on a surface parallel to S can be expressed as a 
function of the distance t between S and this surface as 
r(£) = r ± Ic + 0(£c) 2 , where the sign depends on the 
orientation. We choose the minus sign here, keeping in 
mind that the second monolayer then has the plus sign 
(see figure [T|). Let us now consider the surface such that 
£ = e. If / is written as a function of r(e) = r„ and 
of the curvatures, it features no coupling between these 
variables. This corresponds to the definition of the neutral 



surface |37j , which means that e is the distance between S 
and the neutral surface of the monolayer. On this surface, 
the density which minimizes / for any given membrane 
shape is po (at first order in e). 




Fig. 1. Schematic drawing of a lipid bilayer. The principal 
curvatures ci , C2 and the scaled densities r of both monolayers 
are defined on S. The distances between S and the neutral 
surfaces Af of monolayers ± are denoted e . If the orientation 
convention is chosen in such a way that c < on the drawing, 
the densities on A/' ± are = r ± ±e ± c+C(e 2 ) . In this example, 
monolayer + is constituted of two different lipid species. 



Let us examine the case of a plane monolayer with 
fixed total mass. Its equilibrium density p cq = po (1 + r cq ) 
can be obtained by minimizing its free energy per unit 
mass 

/ = oo + fc + 2o v+o(e 3 ) (2g) 

P Po 2po 



with respect to r. We obtain 



eq 



and thus 



0o (po) 



Po 



Pcq 



2a ' 



O 



Po 



Pcq 



(29) 



(30) 



In particular, (To vanishes when po = Pcq- Note that we 
have assumed that our reference density po was sufficiently 
close to p C q for the second-order expansion of / to be valid 
at p e q- 



3.1.2 Two-component monolayer 

Let us now consider the case of a two-component mono- 
layer. The mass fraction cj) of one lipid species must be 
taken into accout in our monolayer model. In order to 
study small variations of cj) around a reference value <po, 
wc introduce a fourth small variable 



= 0(e) 



(31) 



The expansion of / can now be written as: 



/ = (T + (JiV + y V> 2 + ^{r-ec) - - (c + c ip) c 



+ -c lC2 + 0(e 3 ) , 



(32) 
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where all the coefficients depend on 4>q as well as on po ■ It 
is not necessary to include a term in rip in this expansion. 
As a matter of fact, the conservation of the total mass of 
each lipid species in the monolayer entails that the integral 
of pep is a constant as well as the one of p, so a linear term 
in rip would be redundant with the one in ip. 

The equilibrium density p eq for a flat monolayer with 
a fixed total mass and a uniform lipid composition <p = 4>q 
has the same expression (|29[) as in the case of a single- 
component monolayer, but its value depends on 4>q since 
Co and k do. 

Expanding / in terms of small variables relies on the 
assumption that / is analytical. However, in the case where 
<p is very small, the free energy contains a non-analytic 
part, which reads per unit surface er(l + r)<f>hi<f>, with 
a = k^T po/ '/xi, where /Lti is the mass of one lipid of the 
species with mass fraction cp (see, e.g., the entropy of mix- 
ing in Ref. [3E]). In order to include this term, let us write 
/ as: 



/ 



cr + o-iV> + yV 2 + [cr (1 + r) 01n< 



k 

2 (r 



(cq + CoVO c ' 



-cic 2 + C(e 3 



\ 2 
ec) 

(33) 



where the term between double square brackets must be 
taken into account only if <f> is very small, in which case 
ip stands for <p and not for (cj> — <f>o)/(j>o. In the following, 
the double square brackets will always be used with this 
meaning. 



Using the constraint J A dAp^ (1 + ? ,:t ) = M ± yields 



I A 

ML 



1 ± — / dAc. 
A J A 



(36) 



Let us define the relaxed area of a monolayer as 
the area (defined on S) that it would spontaneously adopt 
in the absence of any constraint: 

A% = M±/p± , (37) 

and its actual area A measured on the neutral surface of 
the monolayer, which verifies: 



We may then rewrite (|3"6"|) as X^p^/k^ 
so that (|35[) becomes 

A ± - A ± 

r ± ere = — — ; . 

A 

Thus, after this partial minimization with respect to r , 
the monolayer free-energy density can be written as 



(38) 
(A±-A-)/A, 



(39) 



f± _ k ( ^ ± 
1 ' 2 \ A 



± \ 2 ± ± 
k cr 



±- 



2 2 -c+^-c? + — ciC2.(40) 



The total free energy F b of the bilayer is obtained by 
integrating / + + /~ over S. Defining the bilayer clastic 
constants and spontaneous curvature as 



3.2 Consistency with the ADE model 

The area-difference elasticity (ADE) model [T71H51[Tg] can 
be deduced from this model by considering a membrane 
made of two one-component monolayers (denoted + and — 
as on figure [T|) with fixed total masses M ± , and by elimi- 
nating the densities by minimization |18| . Choosing p^ = 
p^, which implies <7q = 0, we have 



t± i ± i ± \ 2 i K±c o 

f = — [r ±e c) ± —f-c- 



-cic 2 . (34) 



Since c\ and c 2 are defined on the same surface S of the bi- 
layer, and since a single orientation convention is adopted 
for both monolayers, the curvatures are common to the 
two monolayers. Here e + (resp. e~) is the (positive) dis- 
tance between S and the neutral surface of bilayer + (resp. 
— ). The ± sign in front of ensures that two identical 
monolayers forming a bilayer would share the same value 
of Co- Note that the free energy density / defined in the 
previous section corresponds to /~. 

Let us minimize with respect to the free energy of 
each monolayer with fixed total mass M and total area 
A (defined on 5, like ^ and c\ and c 2 ). To take into ac- 
count the constraints associated with the total masses, we 
introduce two Lagrange multipliers A . We thus minimize 
f ± — X ± Pq(1 + r^) with respect to r*, which gives 



(35) 



we obtain 



(41) 



F b = F b A - 



dAc 2 -K b c b dAc + K b dA Cl c 2 , (42) 

J A J A 

where F b A = [k+(A+ - A+) 2 + k~(A~ - Aq ) 2 ]/{2A). If 
the position of S in the bilayer, which was arbitrary until 
now, is chosen in such a way that 



F b A can be written as 



F b 



= L {A _ Ao f + ^ {AA ^ AAo f 



~ K b 
A ' 4A 

with AA = A~ — A + , AA = Aq - A£, and 



(43) 



(44) 



K 



2k+k- 
k+ + kr 



A 



k+ + k- 



(45) 

Condition (|4"3")l expresses that S is the neutral surface of 
the bilayer [5D]. Thus, e + (resp. e~) is the distance be- 
tween the neutral surface of the monolayer + (resp. — ) 
and the neutral surface of the bilayer. 

The total free energy F b of the bilayer ob- 
tained from our monolayer free-energy density, corresponds 
to the ADE model [TTIITS] . Note that one usually replaces 
the term cx (A — Aq) 2 by a hard constraint A = Aq [TP"] , 
while it is important to keep the area-difference elasticity 
term oc (AA - AA Q ) 2 . 
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4 Projected stress tensor 
4.1 Stress tensor of a monolayer 

In order to obtain the projected stress tensor of our mono- 
layer, we need to determine the projected energy density 
f(f,ip,hi,hij) associated with f(r, ip, c±, C2). We now as- 
sume that the monolayer exhibits only weak deviations 
from the plane (x, y). Then, if hi = 0(e) and hij = 0(e), 
Eqs. ([23" ]) -([2l ]) are satisfied because c\ + c 2 = c = V 2 h + 
C(e 3 ) and ac 2 = det(h tj ) + C(e 4 ). 

The projected free-energy density / and the projected 
mass density p are given by 



/ = / V 1 + (Vhf = / + y {Vhf + G(e 3 ) , (46) 



p = p<Jl + (Vh) 2 =p+^(Vh) 2 + 0(e 3 ), (47) 
Hence, defining the scaled projected density as 

_ = p-po =r+ l (v/l) 2 + , 3^ (48) 

Po 2 
we obtain from Eq. ([55)1 

f = *0 + *ii> + ^ 2 + ^ (Vhf + k - (f - eV 2 h) 2 



2 

+ I (V 2 /i) 2 - ^ (co + coil>) V 2 /i + I det(fty) 



+ [CT(l + f)01n0] +0(e 3 



(49) 



We can now calculate the components of the projected 
stress tensor from (fTo*]) and ([IB"]) , noting that pdf/dp = 
(l + f)df/df. We will restrict ourselves to the first or- 
der in e here, but the tangential components of the stress 
tensor are calculated at second order in Appendix [S] We 
obtain 



(T + f7iV — k{r— eW 2 h) - yV 2 /i 



KC 



S zj = aohj + fee dj (f - eV 2 /i) - -djW 2 h 
+ ^ d ^ + 0(e 2 ). 



(50) 



(51) 



Note that the non-analytic term coming from the entropy 
of a mixture, that we have put between double square 
brackets in (j4"9"j) . does not contribute to the stress tensor. 
This is because this term is proportional to p, so it dis- 
appears when one computes / — pdf/ dp. The expression 
of the stress tensor is thus the same whether 4> is small or 
whether it is close to a finite value 4>o- 

For h = 0, f = and ip ~ 0, ([ST))) gives Sij — <jq 
Hence <jq can be interpreted as the tension of a flat mem- 
brane with uniform density po and uniform lipid compo- 
sition (j)Q. It is consistent with the fact that gq vanishes 
for po = Pcq- Besides, ([30]) may now be interpreted as a 
Hookean law for the tension of a flat membrane with no 
inhomogeneities [39] . 



4.2 Comparison with the stress tensor associated with 
the Helfrich model 

Let us compare our results with those coming from the 
Helfrich model, in which the tension a is a phenomeno- 
logical parameter, namely the Lagrange multiplier imple- 
menting the area constraint. For a monolayer with elastic 
constants \n and ^k, tension a and spontaneous curva- 
ture Co, the Helfrich free energy density is / = ct + kc 2 /4 — 
kcqc/2 + KC1C2/2, so that its projected version reads at 
second order in e: 

/ =ff+ £(Vft) a + |(V a /») 2 -^V^+|det(^). (52) 

The corresponding stress tensor takes the form [32] : 



= (*- ^ 2 h) 6 l3 + ^ hij + 0(e 2 ) , (53) 



KC 



ah, - -S,V 2 /i+0(e 2 



'zj " 3 2 3 



(54) 



Comparing Eqs. (|55 ]) -([5i ]l with Eqs. (f50 ]> — (f5T ]) . we find 
that we may write 



E ij =E* j +0(e 2 ), 



(55) 



Szj = Sfi - ed j( r + ( -± + eat ) d& + 0(e 2 ) , (56) 



if we define 



a = (7 + (i^ - k (f ~ eV 2 /i) + C(e 2 ) 



(57) 



Thus, if the scaled lipid composition ip and the scaled 
density on the monolayer neutral surface f n = r — eV 2 /i + 
C(e 2 ) arc both homogeneous, er is a constant. The stress 
tensor then has the same form in our model as in the Hel- 
frich model. But contrary to the Helfrich tension, our a, 
which may be viewed as a dynamical surface tension, can 
feature inhomogenities. In the inhomogeneous case, new 
terms appear in S z j. Thus, our stress tensor extends the 
one associated with the Helfrich model [32] to the case 
where the lipid density and composition are not homoge- 
neous. 



4.3 Stress tensor in the ADE model 

For a one-component monolayer, the components of the 
stress tensor at first order in e ([SX)]) -([5T ]) are explicitly 
given by 



Exx = o- - k [r - e {h xx + h yy )} 



KC 



r - KC ° h 

■^xy — 2 X V ' 



(58) 
(59) 



S Z x = a h x +kefx--{h xxx + h xyy ), (60) 

where f, = dif and k = k + Ike 2 . The other three compo- 
nents follow from exchanging x and y. 
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Even when the membrane exhibits large-scale defor- 
mations, it is possible to express the stress tensor at a 
given point M in the local tangent frame (X, Y) diagonal- 
izing the curvature tensor. Calling X (resp. Y) the prin- 
cipal direction associated with the principal curvature c\ 
(resp. c 2 ), we have hx = hy = hxy = 0, hxx = ci and 
hyy = c 2 at point M. Hence r = r and c = V 2 /i = c\ + c 2 . 
The components of the projected stress tensor read at first 
order: 

£xx = c - k (r - ec) - ^ c 2 , (61) 
£yy = a -k(r- ec) - ci , (62) 
Sxy = S YX = , (63) 
Szx=-^d x c. (64) 



The tangential stress tensor is thus diagonal. 

Choosing po = p cq and using (|39[) . which comes from 
partial minimization of the monolayer free energies with 
respect to r , we obtain at first order for each monolayer 



£ 



A ± - A?, 



XX 



YY 



A 



± 



c 2 . 



An 



A 



± 



■ Cl 



(65) 
(66) 



Summing the contributions from the two monolayers, we 
obtain the stress tensor of a bilayer in the ADE model, 
still at first order: 



yb 

^xx 



2k' 



YY 
b 

ZX 



= 2k 



A-Ap 
A 

b A-A Q 



A 



K b CQ C 2 



f% Cq C\ 



-k dxc . 



(67) 

(68) 
(69) 



where k b , Aq, K b and Cq are defined in Sec. 

In the Helfrich model, the stress tensor of a bilayer with 
clastic constant n b and spontaneous curvature Cq can be 
written in the principal tangent frame from (|53l) and JS! 
It reads at first order in e: 



s xx = o- - K b c b c 2 , 



^XY 

yli 

^zx 



a - K b c\ ci , 



yH 

^YX 



0, 



-n b dxc . 



(70) 
(71) 
(72) 
(73) 



Thus, the stress tensor in the ADE model has the same 
form as the one in the Helfrich model, with a = 2k b (A — Aq) 
In light of the previous section, this equivalence is not sur- 
prising since in the ADE model, = 0, and f n is homo- 
geneous as shown by (|3"9"j) . 



5 Force density in a monolayer 

5.1 Calculation from the projected stress tensor 



unit area p by taking the divergence of ([50)) and (|5Tj) : 
Pi = dj£ij and p z = dj£ Z j. This force per unit area 
coming from the rest of the monolayer plays an important 
part in a dynamical description of a membrane. Indeed, its 
tangential component is a term of the generalized Navier- 
Stokes equation describing the monolayer, while its nor- 
mal component is involved in the normal force balance 
with the external fluid (see, e.g., Seifert and Langer [2"5]). 
We obtain at first order in e: 



Pi 



-kdi 



(74) 



fc^oV^-^HfeV^+^VV. (75) 

Both of these results give back those of Ref. [3S] in the 
particular case of a bilayer constituted of two identical 
one-component monolayers, if the reference density is po = 
p cq . We have thus justified the expression of these force 
densities from the membrane stress tensor and generalized 
them. 

This force density can also be derived from the general 
expressions (|18[) and ([T9l) . Note that (|T8)) indicates that 
p z = —SF/Sh, which is indeed the force taken into account 
in Eq. (3) of Ref. [25]. Besides, applying (fT^j) to the free- 
energy density (|4"9")l with = shows that, in this specific 
case, ^ = — di(8F/Sf)+0(e 2 y This justifies the "gradient 
of the surface pressure" term — V(5F/5r) used in Eq. (4) 
of Ref. US. 



5.2 Direct covariant calculation 

5.2.1 Definitions and notations 

In this section, we will not restrict ourselves to membranes 
undergoing small deformations around the flat shape. In 
general, a membrane can be considered, in a coarse-grained 
description, as a two-dimensional surface embedded in the 
three-dimensional space. The position of a fluid clement 
in the membrane can be described by a three-dimensional 
vector R(u , u 2 ), where u 1 and u 2 are two parameters la- 
belling each fluid element. Mathematically, these param- 
eters are internal coordinates in the two-dimensional sur- 
face, and physically they correspond to Lagrangian coor- 
dinates. 

We are now going to introduce some basic definitions 
and notations used to describe the shape of a surface in 
/Mffcrcntial geometry 40,34,35 . At each point R of the 
surface, it is possible to define two vectors tangent to the 
surface through 



(76) 



where a £ {1,2}. These two vectors are supposed to be 
linearly independent. Thus, 



Now that we have obtained the stress tensor for our mono- 
layer model, we can calculate the corresponding force per 



ti x t 2 

1*1 X to 



(77) 
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is a unit normal to the surface at point R. The metric 
tensor of the surface can be expressed as 

a a p =t a -tp , (78) 

so that the area element of the surface reads 

dA = y /a'd 2 u, (79) 

where a is the determinant of a a p, and dru = dv}du 2 . The 
inverse metric tensor a a ^ is defined by the relation 



„a/3 ra 

a H ap^ — 6 



(80) 



where <5" is the Kronecker symbol. In the last relation, as 
well as in the following, the Einstein summation conven- 
tion is used. We may now define the contravariant tangent 
vectors as 

r = a Q % . (81) 

A complete description of a surface is given by its metric 
tensor (or first fundamental form) and its curvature tensor 
(or second fundamental form) 



n ■ d a tfj = n ■ d a dnR 



(82) 



The principal curvatures c\ and c 2 of the surface are the 
eigenvalues of b% = a ai b~/p, which enables to express the 
total curvature and the Gaussian curvature from the cur- 
vature tensor: 



c = Ci + c 2 = b% , 
c\C2 = dct bp . 



5.2.2 Force density in a monolayer 



(83) 
(84) 



The surface density of internal forces q in a two-component 
monolayer with free energy F = JdA f can be expressed 
as the functional derivative 



1 



q{u x ,u 2 ) = j= 



SF 



y/E 5R(u\u 2 ) 



(85) 



where, as in the previous sections, p is the total mass den- 
sity of lipids, and <f> the mass fraction of one lipid species 
[341I351I4T] . This expression is a consequence of the prin- 
ciple of virtual work: for a small deformation 5R of the 
membrane at equilibrium, the membrane free-energy vari- 
ation reads 



SF = - JdA q - 5R= - J d 2 u ^ q ■ SR. 



(86) 



For the underlying force balance on each fluid element to 
be valid, the virtual deformation SR must be performed at 
constant total mass dm = pdA = p\fa d 2 u and composi- 
tion in each fluid element. Hence, the functional derivative 
in (|85p must be taken at constant p^fa and pip^fa. 

Let us calculate the force density ([55)1 in a monolayer 
with free-energy density 

f = a + + ^-ip 2 + {a (1 + r) 01n0J + - (r - ec) 2 



- (c + cqtJj) c + -c 



-CiC 2 . 



(87) 



This free-energy density corresponds to (|33[) truncated at 
second order in e. The two constraints on the deformation 
SR, S(py/a) = and S(p<fiy/a) = 0, are equivalent to 



faSr + (1 + r)8y/a = 
y/aSiJj = , 



(89) 



where we have used the fact that 1 + r > 0. To enforce 
these two independent local constraints, we use two local 
Lagrange multipliers, A and p. The principle of virtual 
work (1551) then reads 



8F - J d 2 u {X^Sr + {1 + r)5y/a\ + pVaZS^j} = 5W , 
where 



SW 



dAq-SR. 



(90) 
(91) 



We assume that the topology of the membrane is not af- 
fected by the virtual deformation. The Gauss-Bonnet the- 
orem then ensures that S (JdA C1C2) = 0. The transforma- 
tion of the left-hand side of (JUUJ) can be performed along 
the same lines as in Ref. [12 ■ These calculations, which 
are presented in Appendix [Bj yield: 



(92) 



q ■ t a = — k (1 + r) d a (r - ec) 

01 + a 2 tp - -y^cj d a ip , 

°2 ? k 9 \ k o 
[ cr + a\ip + —ip - -r -kr\c--c° 

+ [kc — k (co + CQip) — 2fce r] C1C2 

+ke{l + r)c 2 + keAr- ^Ac+^yAii, (93) 

where A is a shorthand for the Laplace-Beltrami opera- 
tor (1/ '^/a)d a (a al3 \fadfj). The force density q can be ex- 
pressed from its tangential component (|92[) and normal 
component (|93|) as 



q = {q- t a )t a + (q ■ n)n . 



(94) 



We have thus obtained the general expression of the force 
density in a two-component monolayer with free-energy 
density ([57)1 . This expression gives back the one in Ref. [35] 
in the particular case of a one-component monolayer with 
c = 0. 

Note that, in this section, we have used the free-energy 
density truncated at second order ([8T|) as if it were exact. 
The force density q expressed in ([M)) ~([M )l is the one corre- 
sponding to this model, and it contains second and third- 
order terms. This approach is consistent with the one of 
Rcfs. p^3~TIf3~? ] . However, in the present paper, we have 
constructed the free-energy density as a general expansion 
around a reference state, controlling the order in e of this 
expansion. In our approach, if the free-energy density / is 
kept at second order, the force density can be known only 
at first order. 
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5.3 Comparison between the two results 

In the present paper, except in the previous section 15.21 
we have described membranes in the Mongc gauge, i.e., by 
their height with respect to a reference plane. Such a de- 
scription is very convenient to study the membrane small 
deformations around the flat shape. In the Mongc gauge, 
the position of a fluid clement in the membrane is given by 
R(x,y) = (x, y, h(x,y)), where x and y are Cartesian co- 
ordinates in the reference plane and z = h(x, y) describes 
the height of the membrane with respect to the reference 
plane. Obviously, the coordinates {x, y) of a fluid element 
depend on its position in the membrane (in other words, 
they are Eulerian coordinates). 

In Sec. 15.21 we have found the force density q in a 
monolayer, whatever its shape. In our derivation of q, 
the parameters (u , u ) describing the surface were La- 
grangian coordinates, labelling each fluid element. How- 
ever, as the force density is a physical quantity, it does not 
depend on the parametrization of the surface |34j . Thus, 
the expression we have found for q is valid (for each given 
membrane shape) in the Monge gauge. 

We may now compare the force density obtained in 
Sec. 15.21 with the one obtained from the projected stress 
tensor in Sec. 15.11 For this, we shall write explicitly in 
the Monge gauge the general result obtained in Sec. 15.21 
With u 1 = x and u 2 = y, the tangent vectors read in the 
Monge gauge t\ — (1, 0, d x h) and t 2 = (0, 1, d y h). It is 
then straightforward to find the expression of n, a a p and 
a a/3 in the Monge gauge (see, e.g., Ref. [Mj)- Using these 
explicit expressions, and keeping only first order terms in 
e, (|92p and ([9"3")) can be written as: 

cn = -kdi (f-eV 2 /i) + a 1 d i ip + 0(e 2 ) , (95) 

q z = <j V 2 h - | V 4 ^ + ke V 2 f + + C(e 2 ) , (96) 

where i € {x,y}. We notice that, at this order, (p?5j) is 
identical to ([74")) and (|96[) is identical to ([75)1 . Note that 
q is a force density per actual unit area of the monolayer 
while p is a force density per projected unit area. However, 
this difference is irrelevant at first order. 

We have just shown that the force density obtained 
from the divergence of the projected stress tensor is con- 
sistent with the one calculated directly by using the prin- 
ciple of virtual work in covariant formalism. The projected 
stress tensor thus allows to calculate easily both the nor- 
mal and the tangential components of the force density in 
a membrane in the Monge gauge, without having to resort 
to a covariant formulation. 

6 Applications 

The force density in a membrane with lipid density and 
composition inhomogeneities can be used to understand 
qualitatively and quantitatively the dynamics of a mem- 
brane submitted to a local perturbation. In order to illus- 
trate this, we are going to focus on the local injection of a 
reagent close to a membrane (see Fig. [3J) , which modifies 
locally the properties of the membrane. 




Fig. 2. Local injection of a reagent from a micropipette close 
to a vesicle. The lipids in the external monolayer of the vesicle 
will be affected. The injection is sufficiently local for us to focus 
on a small, nearly-plane zone of the membrane. 



6.1 Forces arising from a modification of composition 

Let us consider initially a one-component flat membrane 
at equilibrium with uniform lipid density. Let us assume 
that, at time t = 0, some lipids in the external mono- 
layer of this membrane (monolayer +, as on Fig. \Q, are 
suddenly chemically modified due to a microinjection of a 
reagent close to the membrane (see Fig. [2]). Then there is a 
fraction <j)(x, y), assumed to be small, of modified lipids in 
this monolayer. The force density in the external mono- 
layer at time t = + , just after the injection, when the 
shape and the density have not changed yet, is given by: 

P + =<TiVHyV^e 2 , (97) 

where e z is a unit vector in the z direction. This force 
density corresponds to Eqs. ([74)1 and ([75)) in the case of a 
flat "+" monolayer with uniform density. Hence, modify- 
ing locally the lipids of a monolayer will gcnerically induce 
a shape instability of the membrane. The internal mono- 
layer is not affected by the chemical modification, so the 
force density remains zero in it. 

Let us take the position of the micropipette injecting 
the reagent as the origin of our (x,y) frame. Then, <f> is 
a decreasing function of the radial coordinate r. Let us 
study the case where <p is a Gaussian: 

<j>(r) = 0o exp ■ (98) 

This can represent the field of modified lipids resulting 
from a diffusion of the reagent in the solution surrounding 
the vesicle before it hits the membrane. It is straightfor- 
ward to calculate the corresponding force density. Its nor- 
mal and radial components, nondimensionalized by their 
maximal values, are plotted as a function of r/R in Fig. [3] 

The constants o\ and So that appear in the force den- 
sity arise from the </>-dependence of the free energy per unit 
area of a monolayer (see Eg. [55)) . Physically, they describe 
the change in the equilibrium density and the spontaneous 
curvature of the membrane due to a generic modification 
of the lipids [2H]- Their signs, and thus, those of the force 
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Fig. 3. Nondimensionalized force density in the external mono- 
layer of a membrane just after a local chemical modifica- 
tion, with Gaussian (/>. Plain line: normal component p z /p™ ax . 
Dashed line: radial component p r /p™ ax . 



force density in monolayers "±" caused by the deforma- 
tion is given by: 



p ± = T keV(W 2 h) - -V 4 /ie 2 



(99) 



This force density corresponds to Eqs. ((74)) and ((75]) in the 
case of one-component monolayers with uniform density. 

Let us consider for instance a Gaussian-shaped defor- 
mation towards the exterior, centered on the origin: 



h(r) = ho exp 



2R 2 



(100) 



The normal and radial components of the corresponding 
force density, nondimensionalized by the absolute value of 
their maxima, are plotted as a function of r/R in Fig. SI 



density components, depend on the nature of the modifi- 
cation. In the case where o\ < and cq < 0, which corre- 
sponds to modified lipids favoring a smaller density and a 
larger curvature (in absolute value), the lipids in front of 
the pipette are submitted to a normal force going towards 
the exterior of the vesicle, which will yield a local defor- 
mation of the membrane in this direction. Meanwhile, a 
radial force drives the lipids of the external monolayer to 
flow in the membrane towards larger values of r, due to 
the fact that the modified lipids favor a smaller density. 

This situation describes well the onset of the shape 
instability studied in Ref. [28], which is induced by mi- 
croinjecting a basic solution close to a giant unilamellar 
vesicle. This instability can yield the formation of a mem- 
brane tubule l22l. 



a, 




Fig. 4. Nondimensionalized force density in monolayer "+" 
just after a local deformation, with Gaussian h. Plain line: 
normal component f»z/|p™ ax |. Dashed line: radial component 
p r /|p™ ax |- In monolayer "— ", p z is identical and p r is opposite. 



6.2 Forces arising from a local deformation at uniform 
density 

Besides the composition, another important effect cap- 
tured by our study is the coupling between the mem- 
brane shape and the density. To shed light onto this ef- 
fect, let us consider a locally deformed membrane with 
uniform density (on the bilayer midsurface) . This can cor- 
respond to a membrane which has deformed very rapidly 
from a flat shape, before the density adjusts to the new 
deformed shape. Indeed, the symmetric density (i.e., the 
sum of the densities in the two monolayers) is not cou- 
pled to the deformation, while the antisymmetric density 
is [55,24,22,28]. Thus, intermonolayer friction is involved 
when the density adjusts to the deformation. The associ- 
ated timescale can be quite large, e.g., a few seconds for 
deformations on length scales of order 20 to 100 fim, so 
there is indeed a lapse when a deformed membrane with 
non-adjusted density exists in the experiments described 
in Ref. [2£]. 

To isolate the effect of the shape and density, we will 
focus on a one-component membrane here (note that the 
external monolayer is a two-component one in Ref. [28]). 
Let us take p = p cq as our reference density in both (iden- 
tical) monolayers. Then, Eq. ((50]) ensures that ctq = 0. The 



Since ke > and k > 0, at small r, the lipids are sub- 
mitted to a normal force going towards the interior of the 
vesicle: this will lead to a relaxation of the deformation. 
This is due to the fact that the membrane considered here 
is symmetric, so its equilibrium shape is flat. Meanwhile, 
a radial force drives the lipids of the external monolayer 
to flow in the membrane, for the density to adjust to the 
shape (see Fig. [5]). Indeed, in a curved membrane at equi- 
librium, the density is uniform on the neutral surface of 
each monolayer, and not on the membrane midsurface. 
The radial forces are opposite in the external and in the in- 
ternal monolayer, because the orientations of these mono- 
layers are opposite while they share the same curvature. 
Hence, for the density to adjust to the deformed shape, 
it is necessary that the lipids in one monolayer slide with 
respect to the ones in the other monolayer, which involves 
intermonolayer friction |28j . 

6.3 Dynamics 

The force densities in the membrane are the basis of a 
hydrodynamic description of the membrane (for a review, 
see Ref. [13]). Our work enables to take into account lipid 
density and composition instabilities in such dynamical 
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Fig. 5. Membrane with Gaussian h, that has just deformed 
from a fiat shape. The density has not adjusted to the new 
deformed shape yet. The lipids are thus at equal distance on 
the midlayer (plain line). The equilibrium density for this de- 
formed shape would correspond to lipids at equal distance on 
the neutral surface of each monolayer (dashed lines). The ar- 
rows indicate the direction of the tangential force density in 
the membrane, consistent with Fig. [4] 



studies. In Rcf. [33], we have used the force densities de- 
rived in the present work to describe the dynamics of a 
membrane curvature instability induced by a local injec- 
tion of a basic solution close to a giant unilamellar vesicle. 

More precisely, we have written down generalized Stokes 
equations describing the balance of the forces per unit area 
acting tangentially in each monolayer. These equations in- 
clude the tangential density of forces Eq. ([73]). the viscous 
force density due to the two-dimensional flow of the lipids, 
the tangential viscous stress exerted by the flow of the sur- 
rounding fluid, and the intermonolayer friction. We have 
also written down the balance of the forces per unit area 
acting normally to the membrane, including the normal 
force density given by Eq. (fT5"j) . and the normal viscous 
stress exerted by the flow of the surrounding fluid. Our 
last fundamental dynamical equation is the conservation 
of mass. 

Thanks to these equations, we have described theoret- 
ically the evolution of the deformation resulting from the 
local chemical modification of the external monolayer by 
the basic solution. This description is developed in detail, 
and compared to experimental results in Ref. [25] , 



7 Conclusion 

We have derived a general formula expressing the pro- 
jected stress tensor in a monolayer as a function of the 
monolayer free-energy density, taking into account inho- 
mogeneities in the lipid density and composition. This 
general formula has been applied to a generic monolayer 
model constructed from basic principles. Our model being 
consistent with the ADE model, we have found in partic- 
ular the stress tensor associated with the ADE model. 

In the Monge gauge, the projected stress tensor pro- 
vides a convenient way of deriving the force density in 
a monolayer, which is the basis of a hydrodynamic de- 
scription of a membrane. The result is consistent with a 
direct calculation of the force density from the principle 
of virtual work in covariant formalism. We have shown 
an example of application to the calculation of force den- 



sity in a locally perturbed membrane. These force density 
constitute the basis of the dynamical study in Ref. [28] . 

Moreover, the stress tensor contains more information 
than the force density, since it provides the actual force 
exerted by a piece of membrane along its edge. Indeed, 
the stress tensor associated with the Hclfrich model [3"Tl 
132] has already been used to study various situations, such 
as the boundary conditions on a membrane with a free 
exposed edge [H] , the adhesion of a fluid membrane [55J 
Hrj] . the surface tension of fluctuating membranes [17,48 
and the force exerted by a fluctuating membrane tubule 
[35] . The stress tensor is also a useful tool in the study 
of membrane- mediated interactions |50U51ll30] . Since the 
stress tensor studied in the present paper generalizes the 
one associated with the Hclfrich model to monolayers with 
density and composition inhomogeneities, it may enable to 
extend such applications. 



A Tangential stress tensor at second order 

For the Helfrich model, the tangential components of the 
stress tensor can be obtained at second order from the 
free energy at second order [3"2"|. On the contrary, here, 
because of the term df /dr, one cannot obtain Sij at order 
e 2 without taking into account in / the terms of order e 3 
which depend on f. 

When such third order terms are included, / becomes 

f = <j a + ntf, + y> 2 + pr (1 + r) <f>ki<f>} 

+ - (r - ec) - - (co + c ip) c+ -c + -cic 2 

H -r 2 tb + —rib 2 H r 3 H -r 2 c -ibrc 

2^ 2^ 3 2 2 v 

+ jrc 2 + ~^rc lC2 + 0\e A ) , (101) 

where C(e 4 ) stands for terms of order e 4 , or terms of 
order e 3 independent of r. Indeed, the latter may be dis- 
carded because they will not contribute to the stress ten- 
sor at order e 2 . The new terms in / may be considered as 
originating from a density-dependence of the constitutive 
constants k, K, K, Co, c"o, <J\ and o-i. 
Using the relations 

/ = fy/l + (Vhf = f + Y ( V/l ) 2 + °V) > ( 102 ) 
c = V 2 /i + C(e 4 ) , (103) 
c lC2 = det(hij) + 0'(e A ) , (104) 

f = r + i ( Vhf + G'(e 4 ) , (105) 
we obtain 

/ = ct + y {Vhf + <j 1 ij + ^^ 2 + la(l + f) 01n 0] 
+ \ (f - eV 2 hf - f (co + co^P) V 2 h + \ (W 2 h) 2 
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det (hi j 
(Vh) 



0-(^ h ) 2 + ^ det (h l3 ) 



2 , £a^a_^^ V 2 fe 
2 Y 2 v 



^ V - ^ V 2 /.) f 2 + | f 3 + 0'^) • (106) 



Eqs. CE5]) and jTOBJ yield 

02 — C 2 2 (TQ + fc 



a + a\ip ■ 



■ V + 



docs not vary during the virtual deformation 51?, wc may 
replace the free-energy density (j87| by 



(111) 



/ = / - 2 ClC2 



in our calculations. The free-energy variation during the 
deformation then reads 

SF = Jd 2 ufS^/a + Jd 2 u^/E5f, (112) 

with 



(k + a[iP) r~(k + 2k') r - + !L-JL (V 2 h)' 
ke - 2 ( c ° + c °^) + kc r + — ip 



\7 2 h 



-— det (%) - (t + g ( c o + °o^0 ^-i 



- (h i3 V 2 h - hidjV 2 ti) + ke (h^f - hidjf) 



KCq 



hid^ + 0(e 3 



(107) 



en + - + [ct (1 + r) (1 + In , 



Sf = 

+ [k(r - ec) + [5-0 In 01 ] or + hSc, 
where we have defined 



Ii = — (c - c - c ?/>) -ek(r- ec) 



Sip 
(113) 

(114) 



The principle of virtual work (|86[) may now be written 



where k = k + 2fce 2 as before. The non-analytic term 
present in / at small <p does not contribute to the stress 
tensor, for the same reason as before. 

In the principal tangent frame, the tangential compo- 
nents of the stress tensor at second order are given by 



XX 



o-o + err ip + °" 2 - ° 2 V 2 - (As + o{il>) r 



SW 



dA \ hSc + [k(r - ec) + [ac61n <p\ - A] Sr 



<ri + ^ ~ + fa (1 + r) (1 + In 0)] - /x 



d 2 u 



/-A(l + r) 



(115) 
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(fc + 2fc') y + ( ke + kc' r + ^ ip ) c 



~2 ( c o + 5 oV0 C2 + ke rci 



k' + £' 



0. 



cic 2 , (108) 
(109) 



Syy can be obtained by exchanging c\ and ci in Sxx- 
This tangential stress tensor thus remains diagonal at sec- 
ond order in the principal tangent frame, like the one as- 
sociated with the Helfrich model 1321. 



where SW is given by (|91l) . The variations 8c and Sy/a only 
come from the variation SR of the shape of the monolayer. 
The coupling between SR and Sr and Sip, which comes 
from the constraints ([58 |) -(|59" ]l . has been accounted for 
by introducing the Lagrange multipliers A and fi. Hence, 
SR, Sr and Sip should now be considered as independent 
variations, and the terms in Sr and Sip must vanish for 
(|115[) to be valid for any virtual deformation, yielding 



A = k(r - ec) + [<7c/>ln0] , 



kc 



(116) 



At = o-i +°2ip- -^c+[cr(l + r)(l + ]n0)] . (117) 



B Covariant calculation of the force density 

In this section, we are going to present the main steps of 
our covariant calculation of the force density in a mono- 
layer, which leads to (j92jl - (j93"]) . This calculation follows 
the same lines as the one in Ref. f42]. 

As the Gaussian curvature term contained in the total 
monolayer free energy 



We thus obtain 



SW = J d 2 u (yfahSc + I 2 6y/a) 



where wc have defined 



h = co + cr\4> + y^ 2 - 2 ( c ° + ^oVO c 



+ -(? + - (ec -r)(r + ec + 2) 
4 2 



(118) 



(119) 



F = dAf = 



d 2 u\fa f 



(110) 



Note that the contribution of the term between double 
square brackets has vanished, as in our calculation of S. 
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Thanks to the relations 

<Va = Va t a ■ St a , (120) 
5c = a afi (d a n) ■ Stp - t a ■ 6 (d a n) , (121) 
Sn = -{n-8t )t , (122) 

SW can be expressed only in terms of 5t a . Performing two 
integrations by parts and using the relations |40j 



t a \/3 — b a pn , 

n\ a = d a n = -bapt 13 



(123) 
(124) 



where g\ a denotes the covariant derivative (associated with 
the metric a a p) with respect to u a of a function g defined 
on the surface 00], we obtain 



SW 



dA 



i a0 h - b ap I 2 ) t 



a a& {d a I 2 )n 



SR. 



(125) 



Identifying (|125|l with (jOTj) for any infinitesimal virtual 
deformation SR, we obtain the sought force density: 



[{a^h - V#I 2 ) ta + a** {d a I 2 ); 



(126) 



Performing the covariant derivative with respect to u" in 
([1261) , usm g (fl23)) - (fT24)l and the relations 00] 



b%\ a = d p c, 



b ap b 



0a = C - ClC 2 ; 



(127) 
(128) 



and taking the scalar product of q with t a (resp. n) finally 
leads to flS2J (resp. @). 
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